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We show how the application of a coupling field connecting the two lower metastable states of a 
A-system can produce a variety of new results on the propagation of a weak electromagnetic pulse. 
In principle the light propagation can be changed from subluminal to superluminal. The negative 
group index results from the regions of anomalous dispersion and gain in susceptibility. 
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A series of experiments have demonstrated both subluminal jj]-^) and superluminal HH propagation of light in a 
dispersive medium. The key to these successful demonstration lies in one's ability to control optical properties of a 
medium by a laser field. Harris et. al. 10] suggested how the electromagnetically induced transparency (EIT) |l(]] can 
be used to obtain group velocities v g fll|| much smaller than the velocity of light in vacuum. Early experiments [[^] 
produced values of group index n g = c/v g in the range 10 2 - 10 3 . Hau et. al. |§ could reduce the group velocity 
to 17 meter/sec in a Bose condensate. This was followed by an experiment in Rb vapor demonstrating reduction of 
group velocity to 90 meter/sec [Q and to 8 meter/sec ||. These experiments were based on the fact that EIT not only 
makes absorption zero at the line center but also leads to a dispersion profile fl(]| , |l2]| with a sharp derivative near the 
line center of the absorption line. In a different development Wang et. al. Q demonstrated superluminal propagation 
following the work of Chiao and coworkers [13-16; see also 17]. Wang et. al. used basically stimulated Raman effect 
with the pumping beam replaced by a bichromatic field. This produces two regions of Raman gain with a region 
in between which has the right anomalous dispersion but with a negligible gain p^Jl9]] . In this communication we 
propose a scheme where by changing a knob - an additional coupling field, one can switch the propagation of light 
from subluminal to superluminal. We basically use the A-system driven by a coherent control field, which has been 
extensively discussed in connection with subluminal propagation We apply, in addition, a field (referred to 

as LL coupling field) on the lower levels of the A-system and demonstrate how the application of the lower level 
coupling field can produce regions in optical response with appropriate dispersion profile. The dispersion can change 
from normal to anomalous depending on the intensity of the LL coupling field [p0|-p2[. In addition, under suitable 
conditions the amplification of the light remains negligibly small. 

We consider the scheme shown in the Fig. 1(a). We consider propagation of light pulse whose central frequency u\ 
is close to the frequency of the atomic transition |1) «-> |3) . We apply a control field on the optical transition |1) <-> |2). 
The transition |2) «-> |3) is generally electric dipole forbidden transition. The states |2) and |3) are metastable states. 
We apply a field of frequency lo 3 on the transition |2) <-> |3). The nature of this field will depend on the level structure. 
It could be a microwave field, say, in case of Na or an infrared field in case of 208 Pb. Moreover, it could be a dc field 
if one is considering transparency with Zeeman sublevels ||. Let 2G = 2d\ 2 .E c /Ti and 2f2 be the Rabi frequencies of 
the control field E c and the LL coupling field, respectively. The state |1) decays to the states |3) and |2) at the rates 
271 and 272. For simplicity we ignore all collisional effects though these could be easily included. What is relevant 
for further consideration is the group velocity v g for the pulse applied on the transition |1) «-> |3). The v g is related 
to the susceptibility x 13 {^i) for transition |1) <-» |3) 

(1) 



l + 2 7 r X ' l3 K) + 2™ 1 ^y^ 

where % i3 (wi) is the real part of x 13 (a;i). We assume that we are working under condition such that Im[x 13 (wi)] = 
x'l ,(wi) ~ 0. The susceptibility x 13 (u>i) will depend strongly on the intensities and the frequencies of the control laser 
and the LL coupling field. We concentrate on the group velocity though actual pulse profiles could be easily simulated 
p3[ . This susceptibility X 13 (wi) is obtained by solving the density matrix equations for the A-system of Fig. 1(a), 
i.e., by calculating the density matrix element P13 to first order in the applied optical field on the transition |1) «-> |3) 
but to all orders in the control field and the LL coupling field. By making a unitary transformation from the density 
matrix p to a via 

Pxi = o 12 e , pi 3 = a 13 e y , p 23 = o~23e , (2) 

we have the relevant density matrix equations 

& u = iGcr 21 + ige~ lAit a 31 - iG*a u - ig*e lAit <r 13 - 2( 7i + 72 )ctii , 
(7 2 2 = iG*<j\2 + i£ler 3 2 — iGa 21 — iSl*a 23 + 272CT11 , 

(7i2 = -[71 + 72 + Ti2 - iA 2 ](7i2 + iGa 22 + ige^ lA ^a 32 - iGa\\ - iWai 3 , (3) 
CT13 = -[71 + 72 + r 13 - z(A 2 + A 3 )](7i 3 + iGo-23 + ige~ lA4t o- 33 - ige^^an - i£la 12 , 
^23 = -(r 2 3 - iA 3 )(7 2 3 + iG*ai 3 + itta 33 - ige~ lAit a 23 - i£lo 22 , 



where T's give collisional dephasings, the detunings Ai, A2, A3, A4 and the coupling constant g are defined by 

d\ 3 .E 

Ai = oji - uj 13 ; A 2 = lu 2 - W12 ; A 3 = lo 3 - lo 23 ; A 4 = Ai - A 2 - A 3 ; g = — — — . (4) 



The susceptibility x 13 can be obtained by considering the steady state solution of (3) to first order in the field on the 
transition |1) <-» |3). For this purpose we write 
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a = a + 9 - e' lA ^ a+ + 9 — e iAit a~ + (5) 

7 7 

The 13— element of a + will yield the susceptibility at the frequency ui\ as can be seen by combining Eqs. (2) and (5) 

N\d 13 \ 2 + 

X 13 (wi) = — a? 3 , (6) 

where N is the density of atoms. In the above equations we have, for simplicity, set ji — 7 2 = 7. The group velocity 
can be obtained by substituting (6) in (1). In presence of the LL coupling field it is difficult to obtain algebraically 
simple expressions for x 13 . However Eqs. (3) can be solved numerically. Doppler broadening can be accounted for by 
using lui — > oji — kv ; W2 — * ^2 — kv ; and by averaging over the Maxwellian distribution for velocities. The velocity 

dependence of cu 3 is insignificant and hence dropped. The parameter S = I B Mc i 1 ) is a measure of Doppler width 

in the Maxwellian distribution p(x) oc exp(— ; x — kv. We show a number of numerical results in Figs. 1 and 
2. We notice from the Fig. 1(c) how the group index n g defined via v g = e/n g , changes from large positive values to 
large negative values and back to positive values as the intensity of the LL coupling field is increased. Thus the LL 
coupling field is like a knob which can be used to change light propagation from subluminal to superluminal. We also 
present the behavior of the corresponding susceptibility for parameters corresponding to superluminal propagation in 
the Fig. 1(b). We see that at Ai = 0, the real part of x 13 exhibits anomalous dispersion whereas the imaginary part 
of X13 is fairly flat and negative and is exactly zero at Ai = 0. The anomalous dispersion along with negative flat 
region in the imaginary part of x 13 is especially attractive for superluminal propagation p4j. In Fig. 1(d) we show the 
behavior of a pulse £{t - L/c) = exp [-(f - L/c) 2 /t 2 ] , £{lo) = exp [-(ui - oJof/T 2 ] , Ft = 2, at the output 

of a medium under condition that group index is negative. The Fig. 1(d) shows that there is no distortion of the 
pulse. For comparison we also show the pulse at the output in the absence of the medium. The advancement of the 
pulse due to medium is seen. The difference in the peak positions is in agreement with negative value of the group 
index. In Fig. 2 we show the results for the group index with and without Doppler averaging. It is known from the 
work of Kash et. al. (J] that Doppler broadening was insignificant in the behavior of the pulse propagation through a 
A-system in presence of a control laser. However the situation changes with the application of the LL coupling field 
at |2) «-> |3) transition (with a wavelength ~ 1.3/ito) particularly in the region where group index is negative. In 
Fig. 2 we also show the results for propagation in a much heavier 208 Pb vapor . This was the system used earlier by 
Kasapi et. al. (!]] to demonstrate subluminal propagation. The application of the LL coupling field can lead to the 
superluminal propagation. The results in this case are not sensitive to Doppler broadening because the Doppler width 
parameter S (= 257) is much smaller than the Rabi frequency G (= 2977) °f the pump. We note that the production 
of superluminal propagation depends very much on the nature of the atomic transitions in the system under study and 
the choice of a large number of parameters such as the powers of the control and coupling fields. From our numerical 
results it is clear that we need a large coupling between |2) and |3). For a magnetic dipole transition between the 
states 1 2) and |3) the requirement of power of the LL coupling field is large and, in principle, this can be met by using 
pulsed fields with a pulse width > [i sec. However if |2) and |3) are chosen to be Zeeman levels, then the available 
dc magnetic field can be utilized to change propagation from subluminal to superluminal. Note that for Rb, a Rabi 
frequency of IOO7 implies a magnetic field of the order of 99.3 Gauss. Another possibility would be to consider an 
effective interaction between |2) and |3) via Raman transition using two other laser fields. The choice of the system 
is quite open and we have essentially shown the "in principle" possibility of light propagation from subluminal to 
superluminal. Thus in conclusion we have demonstrated how the A system can produce a variety of new results if we 
apply an additional LL coupling field. In particular we have demonstrated how the application of the LL coupling 
can produce regions of anomalous dispersion with gain and how this results in superluminal propagation of a weak 
pulse of light. 
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FIG. 1. (a) Schematic diagram of three level A-system. (b) Real and imaginary parts of the susceptibility [Eqs. (6) ; 
X 13 fi.7/N|di3| 2 ] at a probe frequency u>i in the presence of control and LL coupling field Q = 57. (c) Variation of group index 
with the Rabi frequency of LL coupling field. The pulse is taken to have a central frequency on resonance with the transition 
|1) <-> |3). The solid curve of (d) shows light pulse propagating at speed c through 6 cm of vacuum. The dotted curve shows 
same light pulse propagation through the medium of length 6 cm with a time delay — 4.39/isec in presence of a LL coupling 
field with Rabi frequency Q, = 57. The pulse width F is 120KHz. The common parameters of the above three graphs for 87 Rb 
vapor are chosen as density N= 2 * 10 12 atoms/cc, G= IO7, A2 = A3 = 0, Ti2 = Ti3 = T23 = 0, 7 = 3tt * 10 6 rad/sec . 
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FIG. 2. Group index variation with the Rabi frequency of LL coupling field. The curves (a) and (b) are for propagation 
in Rubidium vapor with density N= 2 * 10 atoms/cc. Other parameters are chosen as G= 2OO7, Ai = A2 = A3 = 0, 
Ti2 = Ti3 = T23 = 0, -y — 3 * n * 10 6 rad/sec. For the curve (b) the Doppler width parameter 5 is chosen as 1.33 * 10 9 
rad/sec. Curve (c) shows variation of group index n g with the Rabi frequency of LL coupling field in 208 Pb vapor with density 
N= 2 * 10 14 atoms/cc, G= 2977, Ai = A 2 = A a = 0, r i2 = r 13 = r 23 = 0, 7 = 4.75 * 10 7 rad/sec. 



6 



